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Abstract 

According to the standard models of particle physics and cosmology, there should be a background of 
cosmic neutrinos in the present Universe, similar to the cosmic microwave photon background. The 
weakness of the weak interactions renders this neutrino background undetectable with current tech- 
nology. The cosmic neutrino background can, however, be probed indirectly through its cosmological 
effects on big bang nucleosynthesis (BBN) and the cosmic microwave background (CMB) radiation. 
In this BBN review, focused on neutrinos and, more generally on dark radiation, the BBN constraints 
on the number of "equivalent neutrinos" (dark radiation), on the baryon asymmetry (baryon density), 
and on a possible lepton asymmetry (neutrino degeneracy) are reviewed and updated. The BBN con- 
straints on dark radiation and on the baryon density following from considerations of the primordial 
abundances of deuterium and helium-4 are in excellent agreement with the complementary results 
from the CMB, providing a suggestive, but currently inconclusive, hint of the presence of dark radia- 
tion and, they constrain any lepton asymmetry. For all the cases considered here there is a "lithium 
problem" : the BBN-predicted lithium abundance exceeds the observationally inferred primordial value 
by a factor of ~ 3. 
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0.1 Introduction 

According to the standard models of particle physics and cosmology, neutrinos (known and hypothe- 
sized) are produced, thermalized, and contribute to the total energy density in the early, hot, dense 
Universe, regulating the early Universe expansion rate. Indeed, at the time of big bang nucleosyn- 
thesis (BBN), the contributions to the energy density from baryons, dark matter, and dark energy 
are all subdominant to those from the thermal populations of photons, electrons (e^ pairs), and neu- 
trinos. Since the abundances of the elements formed during the first few minutes of the evolution 
of the Universe depend on the competition between the universal expansion rate and the nuclear 
and weak interaction rates, the very good agreement between the BBN predictions and observations 
(see, e.g.,, ref. [H [21 [31 H] for reviews and further references) depends crucially on the early Universe 
thermalization of neutrinos and places restrictions on the presence of too many (or too few) of them 
or, of too much "dark radiation". At present, BBN and the cosmic microwave background (CMB) 
radiation provide the only probes of the cosmic neutrino background. In addition to their contribution 
to the total energy density, electron neutrinos and antineutrinos play a special role in regulating the 
production of 4 He, the second most abundance element in the Universe. An excess of electron neu- 
trinos over electron antineutrinos (lepton asymmetry; neutrino degeneracy) or, vice versa, will change 
the neutron-to-proton ratio during BBN, modifying, mainly, the BBN-predicted primordial helium 
abundance. 

BBN provides a window on the early evolution of the Universe and a probe of particle physics 
(neutrino physics and more) beyond the standard model (SM). The primordial abundances of the 
elements produced in observationally accessible abundances by BBN (primarily D, 3 He, 4 He, 7 Li) 
depend on three fundamental parameters related to cosmology and particle physics: the baryon abun- 
dance (related to the Universal baryon asymmetry), the expansion rate of the Universe at BBN (a 
probe of dark radiation and the cosmic neutrino background) and, any neutrino degeneracy (lepton 
asymmetry) . 

0.1.1 Baryon Density Parameter 

The most obvious of these parameters is related to the abundance of the reactants, the baryons (nu- 
cleons). Although the very early Universe may have begun symmetric between matter and antimatter 
(riB = rag), long before BBN some yet to be determined mechanism involving the interplay between 
particle physics (violation of the conservation of baryon number, violation of C and CP symmetries) 
and cosmology (out of equilibrium evolution) led to a small but crucial local asymmetry between the 
amount of matter and antimatter in the Universe. After nucleon-antinucleon annihilation, the excess 
(nucleons, by definition) survives {n-Q — ng — > n-Q = njsr) and the number of nucleons in a comoving 
volume is preserved up to the present epoch (and far into the future as well). Since the nuclear reac- 
tion rates depend on the nucleon density, which decreases as the Universe expands, it is convenient to 
normalize the nucleon density to the photon density. After annihilation, the ratio of the nucleon 
number density to the photon number density is unchanged as the Universe expands and cool^J BBN 
depends on the baryon density parameter 7710, defined by 

rno = 10 10 77b = 10 10 (n B /n 7 ). (1) 

2 The number of nucleons in a comoving volume is conserved. Entropy conservation guarantees that, after 
annihilation, the number of photons in a comoving volume is also conserved. 
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The present value (t = to, when the photon (CMB) temperature is To = 2.725 K) of the baryon 
density is often measured by comparing the nucleon mass density to the critical mass density (Qb = 
(Pb/ ' Pcritjo) and, the critical mass density depends on the present value of the Hubble parameter, the 
Hubble constant (H = 100 /ikms -1 Mpc" 1 ) [5], 



n B h 2 = r/io/273.9. 



(2) 



Predicting the baryon asymmetry of the Universe is one of the key challenges confronting the search 
for new physics beyond the standard model. BBN constraints on t?b can help to identify potentially 
successful models of new physics. 

0.1.2 Expansion Rate Parameter 

The scale factor, o = a(t), describes the evolution of the expansion of the Universe. During the early 
evolution of the Universe the expansion rate, as measured by the Hubble parameter, H = (l/a)da/dt, 
is determined by the total energy density which, during those epochs, is dominated by the contributions 
from massless or extremely relativistic particles, "radiation" (R). 



H 2 = 8irGp/3, 



(3) 



where G is Newton's gravitational constant and p = pr. New physics may lead to pr — > p R (dark 
radiation) or, to a modification of the cosmology (general relativity) G — > G' , replacing the SM 
expansion rate with H — )• H' = SH. The expansion rate factor, S, quantifies any departure from the 
standard models of particle physics and/or cosmology. 

Prior to the start of BBN and prior to annihilation (e.g., m e < T <C m^) the only relativistic 
SM particles present are the photons (with g 1 = 2 degrees of freedom or helicities), the pairs 
(g e = 4), and the N^ = 3, left-handed neutrinos and their right-handed antineutrinos (g u = 2N V ), so 
that pr = p 7 + p e + p u . The evolution of the Universe can be scaled out by comparing the total energy 
density to the energy density in the CMB photons. Prior to annihilation, T 7 = T e = T u , so that 
accounting for the different contributions to pr from relativistic fermions and bosons, 
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for = 3. The contribution from possible dark radiation (e.g., sterile neutrinos) may be ex- 
pressed in terms of an equivalent number of SM neutrinos, AN^ [6]. At BBN, which begins prior 
to e^ annihilation, Nj, = 3 + AN^. In this case 
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Allowing for dark radiation, the expansion rate factor, S, is directly related to AN^, 
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It should be kept in mind that new physics (Sbbn 7^ 1) may manifest itself as Gbbn 7^ Go instead 
of AN^ 7^ 0. In this case, comparing Gbbn when T > m e to its present value, 

Gbbn/Gq = Sbb N = 1 + 0.163AN„. (8) 

After annihilation the only relativistic SM particles present are the photons and the neutrinos. 
The SM neutrinos decouple prior to annihilation, when T ~ 2 — 3 MeV, so that when the pairs 
annihilate, the photons are heated relative to the neutrinos. On the assumption that the neutrinos 
are fully decoupled at e ± annihilation, T u /T~, = (4/11) 1 / 3 and, for the SM (AN U = 0), 

However, in the presence of dark radiation or, "equivalent neutrinos" (decoupled, with T = T v 7^ T 7 ), 

s 2 = Br = 1 + (—j—) I (^-) 3 an v = 1 + 0.135AN,,. (10) 



p R V.l-681 y 8 vii 

Since the SM neutrinos aren't fully decoupled at annihilation, they do share some of the energy 
(entropy) when the pairs annihilate [7]. This has the effect of increasing the relative contribution 
of the neutrinos to the total radiation density so that after e^- annihilation, = 3 + AN^ — > N e f j = 
3.046+ AN y . As a result, later in the evolution of the Universe (e. g., at recombination), /)r/p 7 — )• 1.692 
and p' R /p-y 1.692 + 0.227AN,,, so that for T < m e , 

S 2 = p ^ = l+( — j— ) I ( —V 3 AN„ = 1 + 0.134AN,,. (11) 



p R V 1.692/ 8 V 11 

Of course, this post-BBN relation between the expansion rate (S) and the equivalent number of 
neutrinos (AN U ) is only relevant for those epochs when the Universe is radiation dominated. 

BBN codes track the evolution of S from T > m e , prior to annihilation, to T < m e , well after 
e^- annihilation has ended. Since it is important for BBN to follow the evolution of the neutron to 
proton ratio beginning when T > few MeV, prior to annihilation, 

Sbbn = (1 + 7AN^/43) 1//2 = (1 + 0.163AN,,) 1 / 2 . (12) 

A BBN constraint on S is equivalent to one on ANj, (or, on the ratio of Gbbn to its present value, Go) 
and, later in the evolution of the Universe, N e jj = 3.046+AN^. A BBN determination that AN^ differs 
from zero at a significant level of confidence can provide evidence for new physics (dark radiation) 
such as the existence of one, or more, sterile neutrinos (thermally populated) or, a modification of the 
equations describing the expansion rate of the early Universe (S'bbn 7^ !)■ 



0.1.3 Neutrino Degeneracy Parameter 

Since the charge neutrality of the Universe ensures that any electron excess is tied to the proton 
excess (the baryon asymmetry) , a non-zero lepton asymmetry much larger than the baryon asymmetry 
(t/b ^ 10 -9 ) must be hidden in the neutrino sector. An excess of neutrinos over antineutrinos (or, vice- 
versa) requires a non-zero neutrino chemical potential, p u . The dimensionless degeneracy parameter 
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is the ratio of the neutrino chemical potential to the neutrino temperature, = fi u /T u ; £ u is preserved 
as the Universe expands and cools. In analogy with the parameterization of the baryon asymmetry 
by ?7b = ( n B — n g)/ n 7 ~~ ^ n B/ n 7) a lepton (neutrino) asymmetry may be parameterized by 



7T 
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(13) 



where the sum is over the three SM neutrino flavors (a = e, fj,, r). Generally, mixing among the SM 
neutrinos ensures that the three chemical potentials are equilibrated. In the following it is assumed 
that £ = £ e = £ M = £ r . In this case, 
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(14) 



An asymmetry between electron neutrinos and electron antineutrinos has a direct effect on BBN 
through the charged current weak interactions which regulate the neutron-to-proton ratio (p + <B- 
n + v e , n + e + p + z/ e , n p + e~ + D e ) (see, e.,g., [U El HOj HU H21 [131 E! anci further references 
therein) . Since the relic abundance of 4 He depends directly on the neutron-to-proton ratio when BBN 
begins (and during BBN), it provides a sensitive probe of any lepton asymmetry. The abundances of 
the other light nuclides produced during BBN are less sensitive to £. 

A subdominant effect (usually) of a non- negligible neutrino degeneracy (t/l 3> Vb) is to enhance 
to the contribution of the neutrinos to the early Universe energy density. This is equivalent to a 
contribution to AN^ where, for £ e = ^ = £ T = £, 
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(15) 



Note that for |£| < 0.1, AN ;/ (£) < 0.013, which is likely small compared with anticipated uncertainties 
in ANj, inferred from BBN or the CMB. 

At present and, likely for the foreseeable future, BBN provides the only window to a universal 
lepton asymmetry. 



0.2 The BBN Predicted Abundances 

For BBN within the context of the standard models of particle physics and cosmology (SBBN), along 
with some well defined extensions of them, only the light elements D, 3 He, 4 He, and 7 Li are produced 
in observationally interesting abundances. The BBN-predicted relic abundances of these light nuclides 



depend on the three fundamental parameters introduced in §0.1 [101 H21 [TBI EE CI EH E2]- Over 
limited, but interesting ranges of these parameters, the results for the abundances of these nuclides 
extracted from numerical BBN codes, are well fit (within the quoted errors) by |13 [ \T\ IT4"1 115j . 

y D p = 10 5 (D/H) P = 2.60(1 ± 0.06)(6/r/ D ) 16 = 45.7(1 ± O.Oe)^ 16 , (16) 

Y P = 0.2477 ± 0.0006 + 0.0016(?7hc - 6) = 0.2381 ± 0.0006 + 0.00167^, (17) 
y LiP = 10 10 (Li/H) P = 4.82(1 ± 0.10)(77 Li /6) 2 , A(Li) = 12 + log(Li/H) P , (18) 
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where 

VD = mo ~ 6(5 - 1) + 5£/4, (19) 

me = mo + 100(5 - 1) - 575£/4, (20) 

m = mo - 3(5 - 1) - 7£/4. (21) 



The relation of r/He to £ in Eq. 20 is the one that appears in Kneller & Steigman (2004) [13] . An 
inadvertant typo in Steigman (2007) [TJ was propagated in Simha & Steigman (2008b) [15]. The very 
small difference this typo generated, 575 vs. 574, has no effect on the quantitative results presented 
in those papers. In the above equations, (D/H)p and (Li/H)p are the ratios by number of deuterium 
and of lithium ( 7 Li) to hydrogen respectively, and Yp is the 4 He mass fraction. 



There are some small but interesting changes in the numerical values in Equations 16-18 from 



earlier versions of these relations [131 El EH ES] j Equations 19-21 are unchanged. In a recent paper, 
Nollett & Holder (2011) [4] called attention to the tension between the experimental and theoretical 
determinations of the d(p, 7) 3 He cross section, important for predicting the BBN deuterium abundance. 
Nollett & Holder argue for preferring the theoretical calculation over the experimental result which, 
they suggest, may be affected by a normalization error. Adoption of the theoretical calculation results 
in a ~ 6% reduction in the BBN-predicted D abundance. I have preferred to "split the difference", 
reducing the previously predicted abundance by 3% but, doubling the error uncertainty (the error in 
the BBN-predicted value of D/H for a fixed value of t?d) from 3% to 6%, resulting in the numerical 



values shown in Eq. 16 



Given the role of the neutron-to-proton ratio at BBN on the predicted relic abundance of 4 He, Yp 
depends, albeit weakly, on the value of the neutron lifetime (mean life). Quite recently, the Particle 
Data Group [17J, in response to discrepant experimental data, decided to change its recommended 
value for the neutron lifetime from r n = 885.7 ± 0.8 s to r n = 881.5 ± 1.5 s. This change results in a 
small but noticeable reduction in the predicted value of Yp by 0.0008 and, a very small increase in the 



associated uncertainty in the predicted value of Yp, 0.0005 — > 0.0006. Also incorporated into Eq. 17 
is the Mangano et al [7] correction to the helium abundance resulting from the incomplete decoupling 
of the neutrinos at annihilation. 

The intense interest in recent years in the "lithium problem(s)" has led to an extensive reevalu- 
ation of the relevant nuclear reaction rates |181 [T9l [20] , leading to an increase in the BBN-predicted 
abundance by ~ 12% from the result presented in Steigman 2007 [lj, further exacerbating the lithium 
problem to be discussed below. This is reflected in Eq. 18. The d(p, 7) 3 He rate which plays a role 



in the primordial deuterium abundance also impacts the BBN-predicted lithium abundance [I] (K. 
Nollett, Private Communication). This is taken into account in our error estimate. 

It should be emphasized that these fits are not analytic approximations to the results from a 
numerical BBN code. Rather, they are fits to the results from such a code which are primarily 
simple, involving a minimal number of numerical values, chosen to only one or two significant figures. 
For limited but interesting ranges of the parameters (5.5 < 7710 < 6.5, 0.85 < 5 < 1.15 (—1-7 < AN,, < 
2.0), —0.1 < £ < 0.1) these fits agree with the numerical results from this and other codes within the 
quoted uncertainties. 

It is worth noting that the BBN-predicted abundances of D and 7 Li are mainly sensitive to the 
baryon abundance while that of 4 He is more sensitive to non-standard physics (AN„ and/or £). How- 
ever, both D and 7 Li are weakly dependent on ANj, and/or £ and, 4 He is weakly dependent on r/p. 
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0.3 The Observationally Inferred Primordial Abundances 

Any conclusions about new physics (e.g., AN„ / 0?, £ 7^ 0?) based on BBN depend on the adopted 
primordial abundances. The relic abundances of D and 4 He are key to the conclusions reviewed and 
updated here. It may have been noticed that 3 He failed to be included in the discussion in the previous 
section. The reason is that 3 He is only observed in the interstellar medium of our galaxy, which consists 
of gas that has been processed through many generations of stars. The large and uncertain corrections 
for post-BBN stellar processing make it difficult to infer the primordial abundance of 3 He using the 
current data (see, e.g., [21J for discussion and further references). In addition, the observationally 
inferred abundances of lithium in the oldest, most metal poor stars in the Galaxy are systematically 
lower, by factors of ~ 3 — 4, than the BBN-predicted values (the lithium problem or, one of several 
lithium problems). Whether this discrepancy results from poorly understood corrections for stellar 
structure and/or evolution or, is a hint of new physics, remains unclear at present. In the confrontation 
of the BBN predictions with the observational data only D and 4 He will be used to constrain various 
combinations of rjxo, AN^, £, and the results will be used to predict the primordial lithium abundance 
which will then be compared to the observations. For a much more detailed, albeit not entirely up 
to date, discussion of the observational data and, in particular, the problems (real and potential) 
associated with them, see my recent reviews [HE]. 

0.3.1 Primordial Deuterium 

In the post-BBN Universe, as gas is cycled through successive generations of stars, deuterium is de- 
stroyed, not produced [22]. The post-BBN evolution of deuterium is simple and monotonic (decreasing 
abundance). As a result, observations of deuterium at high redshifts (z) and/or of gas at low metal- 
licity (Z), where very little of the primordial gas has been cycled through stars, should provide a view 
of very nearly primordial deuterium. While interesting on their own, observations of deuterium in 
the chemically evolved Galaxy or the solar system are of relatively little use in constraining the relic 
deuterium abundance. Observations of D at high z and low Z are provided by the QSO absorption-line 
systems (QSOALS) [231 EU E51 [26] [27] . Since deuterium is observed by absorption of background QSO 
light in the wings of the much larger hydrogen absorption, exquisite velocity information about the 
absorbing gas is crucial to a meaningful determination of the D/H ratio. This, and other contributors 
to potential systematic errors, has limited the number of "robust" D abundance determinations from 
high— z, low— Z, QSOALS. In Fig.[T]are shown 12 high— z, low—Z D/H determinations as a function of 
the absorbing redshift (upper left panel), of the metallicity (upper right panel), and of the Hi column 
density (lower panel). The open symbols reflect a subjective judgment of abundance determinations 
which may be more uncertain than indicated by their error bars (perhaps all the data should be plotted 
with open symbols). 

For all 12 D abundance determinations, the weighted mean abundance is < logyo > = 0.42. How- 
ever, it is clear (by eye) that there is an excessively large dispersion among the individual abundance 
determinations (e.g., the reduced \ 2 f° r H degrees of freedom is x 2 /dof = 4). It is also clear from the 
three panels in Fig. [T] that the spread in abundances does not correlate with either redshift, metallicity 
or, H 1 column density. The absence of any correlations suggests that it is unlikely that the spread in 
the observed deuterium abundances results from post-BBN evolution. Lacking a well motivated under- 
standing of the cause(s) of the observed dispersion, various statistical approaches have been adopted 
for estimating the uncertainty when identifying the weighted mean abundance from the observations 
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Figure 1: The log of the deuterium abundances (^d = 10 5 (D/H)), and their 1 a uncertainties, inferred 
from 12 low— Z, high— z QSOALS. In the upper left panel t/d is shown as a function of the redshift. In 
the upper right panel is shown as a function of the metallicity (squares for silicon and triangles for 
oxygen). In the lower panel j/d is shown as a function of the neutral hydrogen column density. The 
filled symbols reflect a subjective determination of the more robust determinations of D/H compared 
to the open symbols. The blue solar symbol in the upper right panel is the solar system (pre-solar 
nebula) D abundance 
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with the primordial deuterium abundance. Consistent with these more sophisticated approaches, if 
the quoted errors for each of the data points is simply doubled (lowering the reduced x 2 /dof from 4 
to 1), it leads to the following estimate of the relic abundance, 

log y DP = 0.42 ±0.02 (y DP = 10 5 (D/H) P = 2.63 ± 0.12). (22) 

This result is consistent with that quoted in Pettini & Cooke (2012) [27], who used a slightly different 
set of D/H observation^} Adopting this estimate for the primordial D abundance, Eq. 16 results in 



7/d = 5.96 ± 0.28. (23) 

0.3.2 Primordial Helium 
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Figure 2: Helium abundance (mass fractions, Y) determinations from the sample of extragalactic 
Hn regions studied by Izotov & Thuan (2010) |29j as a function of the corresponding oxygen abun- 
dances (O/H by number). The solid line is the Izotov & Thuan best fit to a linear Y versus O/H 
correlation. 



3 For their new, most precise individual deuterium abundance determination, Pettini & Cooke [27! find yr> = 2.53±0.05 
which, if identified with the primordial deuterium abundance, corresponds to r/u = 6.10 ± 0.24. 
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As was the case for deuterium, the post-BBN evolution of helium ( He) is simple and monotonic. 
As gas is cycled through stars, hydrogen is burned to helium (and beyond) and the helium abundance 
increases with time and with metallicity. The strategy, therefore, is to concentrate on determining the 
helium abundance in the most nearly primordial regions of low metallicity which, as with deuterium, 
lie outside of the Galaxy. In Hn regions, regions of hot, ionized gas, recombinations of hydrogen and 
helium result in observable emission lines which can be used, along with models of the H n regions and 
a knowledge of the associated atomic physics, to infer the helium abundance. The current inventory of 
helium abundance determinations from relatively low metallicity, extragalactic H n regions approaches 
~ 100 (Izotov & Thuan 2010 (IT) [29]). These data for the inferred helium mass fraction (Y) as a 
function of the corresponding oxygen abundance (O/H by number) are shown in Fig.[2j The sheer 
size of this data set leads to relatively small, formal statistical errors, magnifying the importance of 
taking proper account of the many possible sources of systematic errors. While some have employed 
a posteriori selected subsets of the IT data for more detailed analyses, the sources and magnitudes of 
systematic errors have rarely been addressed. As a result, the uncertainty in the inferred primordial 
helium mass fraction is currently dominated by systematic errors (both the known unknowns and the 
unknown unknowns). As a result, estimating the size of the true uncertainty in the observational 
determination of Yp is largely guesswork. From a linear fit to their Y - O/H data, IT find the 
intercept, providing an estimate of the primordial helium abundance, 

Y P = 0.2565 ± 0.0010 (stat) ± 0.0050 (syst). (24) 

In the analysis presented here the statistical and systematic errors are, arbitrarily, combined lin- 
early, leading to the estimate of the primordial abundance adopted here, 

Y P = 0.2565 ± 0.0060. (25) 

This estimate of the relic 4 He abundance is consistent with other, recent estimates based on analyses 
involving limited subsets of the IT data. In some of those other analyses, a linear Y versus O/H fit is 
forced on data which is consistent with no correlation between Y and O/H. Not surprisingly, the result 
of such analyses is a slope which is consistent with zero at less than 1 a. However, the large uncertainty 
in the slope inferred from such fits leads to an estimate of the intercept (Yp) with excessively large 
errors, which have nothing to do with either the statistical or systematic errors. The errors simply 
reflect the uncertainty in the slope for uncorrelated data. Even worse, since these fits are consistent 
with an unphysical, negative Y - O/H slope at < 1 er, they lead to an unphysical upper bound to Yp. 



In combination with Eq.[T7j the relic abundance adopted here results in, 

r/He = 11-50 ±3.77. (26) 



0.3.3 Primordial Lithium 

Compared to the post-BBN evolution of D and 4 He, the evolution of lithium is more complicated 
and uncertain, similar to that of 3 He. As gas is cycled through stars, most of the pre-stellar lithium 
is destroyed. However, some lithium may avoid nuclear burning if it remains in the cooler, surface 
regions of the coolest, lowest mass stars. Observations suggest that some stars (the "super-lithium 
rich" red giants), during some part of their evolution, are net producers of lithium, although it is not 
entirely clear if such stellar produced lithium is returned to the interstellar gas before being destroyed. 
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Finally, it is well known that collisions in the interstellar medium between cosmic rays, mainly alpha 
particles, and interstellar gas nuclei, primarily CNO nuclei, break up those nuclei producing lithium 
( 7 Li, along with 6 Li and, isotopes of Be and B). The net effect of post-BBN production, destruction, 
and survival is difficult to model precisely. However, there is observational evidence supporting a 
lithium abundance which increases along with the heavy element abundance (metallicity) , suggesting 
an overall increase of the lithium abundance with time. 

As we are interested in samples of the most nearly primordial material, the best (only) targets 
for determining the relic lithium abundance are the oldest, most metal poor stars in the Galaxy. 
If the metallicity is sufficiently small, so that the material in these stars has suffered very little 
processing, observations should find a "lithium plateau". That is, for such metal poor stars the lithium 
abundance should be uncorrelated with metallicity, revealing the primordial lithium abundance. Since 
the observationally inferred lithium abundance (relative to hydrogen) is so small (by number, Li/H 
~ 10~ 10 — 10~ 9 ), it is common to measure it on a logarithmic scale by the quantity A(Li) = 12 + 
log(Li/H). The metallicity is usually quantified by comparing, also on a logarithmic scale, the iron 
abundance (Fe/H) to that in the Sun: [Fe/H] = log(Fe/H) — log(Fe/H) . Observations of stars with 
—2.5 < [Fe/H] < —1.0 (those with ~ 0.3% to ~ 1% of the metallicity of the Sun) do appear to lie 
on a plateau, the "Spite plateau" [30], at a level of A(Li) 2.2 ± 0.1 As will be seen below, all 
of the BBN predictions are close to A(Li) ~ 2.7 or, a factor of ~ 3 higher than suggested by the 
observations [32l I33| 13^] . This is one of the lithium problems. However, as recent observations of even 
more metal poor stars (—3.5 < [Fe/H] < —2.5) have accumulated, the lithium plateau appears to 
be transforming into a "lithium cliff", with lower lithium abundances (A(Li) ~ 2.1 ±0.1) correlating 
with lower metallicities j.32[ [33], |3~^] . This trend is puzzling and not understood at present; another 
lithium problem. Since this second lithium problem is currently unresolved, it is not clear if the 
lower abundances are too be interpreted as suggesting an even lower value for the primordial lithium 
abundance, further exacerbating the original lithium problem or, if they are telling us something about 
the evolution of the oldest, most metal poor stars in the Galaxy which will require us to reevaluate both 
lithium problems, along with the observationally inferred value of the primordial lithium abundance. 

For these reasons (the lithium problems), lithium does not provide a useful probe of BBN at present. 
For a recent review of the current lithium data and possible resolution of the lithium problem(s), see 
Fields 2011 [35]. In the following, BBN, possibly in combination with the CMB, will constrain the key 
parameters using the observationally inferred primordial D and 4 He abundances and, those parameter 
combinations will be used to predict the BBN abundance of lithium (reinforcing the problem of the 
low observed abundances). 

0.4 BBN Constraints On The Fundamental Parameters 

The BBN-predicted primordial abundances depend on all three parameters: 7710, ANj,,£. Because of 
the uncertainty of if, or how, to use lithium ( 7 Li), here we limit ourselves to employing only the D and 
4 He abundances. Then, without recourse to additional, non-BBN data, only two of these parameters 
can be constrained at a time. Here, we will consider constraints on the baryon density (r/10) from BBN 

4 Recently, Nissen & Schuster [31] have suggested that stars with —1.5 < [Fe/H] < —0.7 were formed with a lithium 
abundance close to primordial, which they estimate as A(Li) = 2.58 ± 0.04 (stat) ± 0.10 (syst), much closer to the BBN- 
predicted abundances discussed below. They attribute the lower observed abundances to lithium depletion in the stellar 
atmospheres. 
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(using D as our primary baryometer) and on the parameter pairs {7710, AN,} and {7710, £} (using D 
and 4 He in combination) and, we will comment on the result of using complementary data on 7710 or 
AN, from the CMB in order to constrain all three parameters simultaneously. 



0.4.1 BBN Constraint On The Standard Model Baryon Density (AN, = = 

Before entertaining the possibility of new physics (dark radiation and/or lepton asymmetry), the D 
abundance may be used to provide a standard model BBN (SBBN) constraint on the baryon density; 
the SBBN-predicted D abundance is much more sensitive to the baryon density than is the helium 
abundance. Assuming that AN, = = £, rj D = r] 10 = 5.96 ± 0.28 (n B h 2 = 0.0218 ± 0.0010), in 
excellent agreement with the value found from the CMB (e.g., from WMAP |36j). 

This SBBN-inferred value of the baryon density may be used to predict the primordial abundances 
of the other light nuclides, 3 He, 4 He and 7 Li. There is good agreement between the SBBN-predicted 
value of the 3 He abundance and an upper bound to it inferred from observations stellar-processed gas 
in the Galaxy [21] ■ For SBBN the relic 4 He mass fraction is predicted to be Y P = 0.2476 ± 0.0007, 
a value smaller than the observationally-inferred abundance adopted here but, within ~ 1.5 ex of it 



(see §0.3.2). The problem for SBBN is lithium. For the above value of the baryon density and for 
AN, = = £, the predicted primordial lithium abundance is A(Li) = 2.68 ± 0.06. This is higher, by 
a factor of ~ 3 — 4, than the values inferred from observations of the most metal poor (most nearly 
primordial) stars in the Galaxy. This is (one of) the lithium problem(s). 

0.4.2 BBN Constraints On The Baryon Density (7710) And Dark Radiation (AN,) 



If it is assumed that there is no lepton asymmetry, £ = 0, Eqs.[19]& 20 may be solved for 7710 and S 
(AN,) in terms of rjn and r/He- 

106(5 - 1) = r? He - Vd = 5.54 ± 3.78, (27) 

106r/io = IOOt/d + 6r/H e = 665 ± 36. (28) 

In Fig. [3] are shown contours of constant values of 7710 and AN, in the Yp — yop plane, along with 
the observationally inferred values of ?/dp an d Yp and their 1 a error bars. Measurements of ypjp an d 
Yp constrain 7710 and AN,. From BBN, using D and 4 He, it is found that 7710 = 6.27 ± 0.34 (Q^h 2 = 
0.0229 ± 0.0012) and AN, = 0.66l[j; 4 £ (N e// = 3.71±° H). Figfl shows the 68% and 95% confidence 
contours corresponding to these results. While AN, ~ 1 (a sterile neutrino?) is somewhat favored, 
this result is also consistent with no dark radiation (AN, = 0) within 95% confidence. However, as 
may be seen from Fig.|4j the presence of two sterile neutrinos is disfavored at > 95% confidence. 

The observationally-inferred values of 770 and 7/He may be used to predict the primordial abundance 
of lithium synthesized during BBN. 

10677Li = 103t? d + 37/He = 648 ± 31. (29) 

This leads to the prediction that A(Li) = 2.70 ± 0.06, far in excess of the lithium abundances inferred 
from the observations of the most metal-poor stars in the Galaxy, reinforcing the SBBN result of a 



lithium problem (see §0.4.1). 
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10 5 (D/H) p 

Figure 3: Contours of constant values of AN^ (red) and 7710 (blue) in the Yp — yop plane. From bottom 
to top the red curves correspond to AN y = 0, 1,2. From left to right the blue curves correspond to 
7710 = 7.0, 6.5, 6.0, 5.5. Also shown (filled circle and error bars) are the adopted primordial abundances 
of D and 4 He and their 1 a uncertainties. 

0.4.3 BBN Constraints On The Baryon Density (r/i ) And Lepton Asymmetry (£) 

If it is assumed that there is no dark radiation (ANj, = 0), the observationally-inferred abundances of 
D and 4 He may be used to constrain the baryon density and any lepton asymmetry. 

145£ = 77d - r/He = -5.54 ± 3.78, (30) 

1167710 = H57/D +77 He = 697 ±32. (31) 

In Fig. [5] are shown contours of constant values of r/io and £ in the Yp — j/dp plane, along with the 
adopted values of yop and Yp and their 1 a error bars. Measurements of yrjp and Yp constrain 
7710 and £. From BBN using the adopted primordial D and 4 He abundances it is found in this case 
that 7710 = 6.01 ± 0.28 {n B h 2 = 0.0219 ± 0.0010) and £ = -0.038 ± 0.026. The latter result is 
consistent with £ = at ~ 1.5(7. At 2 a, this result provides an upper bound to the magnitude of 
the neutrino degeneracy parameter (|£| < 0.090) which can be used to constrain the contribution to 
ANj, resulting from the presence of the "extra" energy density associated with an excess of neutrinos 
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Fi gure 4: The BBN-inferred 68% (solid) and 95% (dashed) contours in the AN;, — 7710 plane derived 
from D and 4 He assuming that £ = 0. 



over antineutrinos or, vice-versa: ANj,(£) < 0.011 at ~ 2a. Fig. [6] shows the 68% and 95% confidence 
contours corresponding to these results. This result is consistent with no lepton asymmetry to better 
than 95% confidence. 



As before (see § 0.4.2 ) the observationally-inferred values of t/d and ?7He may be used to predict the 
BBN lithium abundance 



145r/Li = 1427? D + 37/hc = 881 ± 41. 



(32) 



This leads to the prediction that A (Li) = 2. 69 ±0.06, in almost exact agreement with the corresponding 
prediction in the presence of dark radiation (assuming that £ = 0), reinforcing, once again, the lithium 
problem. 



0.5 Discussion 

The simplest (least interesting?) assumption is that of no new physics; the standard model with 
no dark radiation (e.g., no sterile neutrinos) or a significant lepton asymmetry (ijb <C ?]l "C 1). In 
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Figure 5: Contours of constant values of £ (red) and t]iq (blue) in the Yp — t/dp plane. From bottom to 
top the red curves correspond to AN,, = 0, —0.05, —0.10. From left to right the blue curves correspond 
to r/io = 6.5,6.0,5.5. Also shown (filled circle and error bars) are the adopted primordial abundances 
of D and 4 He and their 1 a uncertainties. 



this case the SBBN-predicted baryon abundance (see § |0.4.l| ), tt B h 2 = 0.0218 ± 0.0010, is in excellent 
agreement with the more precise value found from the CMB, ft B h 2 = 0.0226 ± 0.0004 [36]. For this 
value of the baryon density (and AN^ = £ = 0) the SBBN-predicted abundance of 3 He is consistent 
with the primordial value inferred from Galactic observations and, the relic abundance of 4 He agrees 
with the observationally-inferred value adopted here, within ~ 1.5 a. As is by now well established, 
the predicted relic lithium abundance exceeds the values of the lithium abundance inferred from 
observations of the most metal-poor stars in the Galaxy by factor of ~ 3 — 4. 

Setting aside for the moment the possibility of a large neutrino degeneracy 3> rj B ) but, allowing 
for the presence of dark radiation (ANj, ^ 0), BBN along with the adopted primordial abundances 
of D and 4 He may be used to constrain ANj, and £l B h 2 



-0.47 /at 
-0.45 ( iN e// 



3.71™) and n B h 2 



In this case (see §0.4.2) it is found that 
0.0229 ± 0.0012, in excellent agreement with the 



ANy = 0.66; 

values of these parameters inferred from various CMB observations [36} \37\ [38| 139] . When derived 
from the CMB, the errors on ^ e ff = 3.046 + AN V are larger (typically by a factor of ~ 1.5 — 2), and 
those on Q B h 2 smaller (typically by a factor of ~ 3), than the corresponding BBN uncertainties. This 
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Figure 6: The BBN-inferred 68% (solid) and 95% (dashed) contours in the £ — r]io plane derived from 
D and 4 He assuming that AN„ = 0. 



is likely to change when the PLANCK collaboration analyzes its CMB data. The PLANCK constraint 
on AN„ is expected to be more precise compared to the BBN value by a factor of ~ 2.5, while that on 
the baryon density parameter should be more precise than the BBN value by an order of magnitude 
|414 l42j. The BBN-predicted relic lithium abundance when ANj, ^ is hardly changed from the SBBN 
case, reinforcing the lithium problem. As may be seen from Fig.|4j while the result for AN V (N e jj) is 
closer to AN^ = 1 than to AN y = 0, offering some support for the existence of one sterile neutrino, it 
is consistent with ANj, = at ~ 1.5 a. In contrast, the existence of two sterile neutrinos is disfavored 
by the BBN data at > 95% confidence. 

In Fig. [7] the BBN constraints on AN;, and VL^h 2 are compared with those from recent CMB 
analyses. Here, too, it appears that current data have a preference for one sterile neutrino while being 
slightly inconsistent with two sterile neutrinos {e.g., BBN and SPT in Fig. [7]) or with no dark radiation 
{e.g., WMAP7 and ACT in Fig.0). 

Unlike the CMB, BBN has the potential to probe a non-zero (albeit relatively large) lepton asym- 
metry (neutrino degeneracy). Current data, driven by the adopted 4 He abundance, are consistent with 
a small, negative value for the neutrino degeneracy parameter, £ = —0.038 db 0.026 which, however, 
is only ~ 1.5 a from zero. A more precise result will only come when (if) there is a reduction in the 
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Figure 7: Comparing the BBN predictions of N e jy and Qp^ 2 with those from various CMB deter- 
minations: BBN D + 4 He (red filled triangle), BBN D + WMAP7 [36] n B h 2 (red open triangle), 
WMAP7 [36] (blue filled square), ACT [37j (green filled pentagon), SPT [38] (purple filled circle), 
SPT + Clusters [39J (purple open circle). 



error of the observationally- inferred value of Yp, whose uncertainty is dominated by systematics. 

0.5.1 Sensitivity Of AN^ And £ To Primordial Helium And Its Uncertainty 

The BBN-predicted helium abundance is sensitive to the early Universe expansion rate (S or, dark 
radiation AN„) and to a lepton asymmetry (£) and very insensitive to the baryon density (r/e). 
The results which have been presented here for AN^ and £ are mainly driven by the adopted value 
for the primordial helium abundance, and its uncertainty. But, the observationally inferred helium 
abundance, Yp, is a quantity which has changed dramatically over time as more and better data have 
been acquired and more careful analyses of the data have been performed. In the left panel of Fig. [8] 
is shown a chronology, over the past ~ 20 years, of the published observational determinations of the 
primordial helium mass fraction, revealing a nearly monotonic increase of Yp with time. In the right 
panel of Fig. [8] the chronology of the corresponding ANj, values is shown, mirroring the increase in Yp. 
Notice that only very recently, within the past 5-7 years, do the data begin to favor AN^ > 0. 



z 4 
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Figure 8: The left panel shows a history of the primordial helium mass fraction (Yp) determinations 
as a function of time. The same symbols/colors correspond to determinations from collaborations 
involving many of the same participants and/or the same observational data. The right panel shows 
the corresponding chronology of BBN-determined values of AN,. The dashed line shows the SM 
result, AN, = 0. 



0.5.2 Constraints On AN, From BBN D And The CMB-Inferred Baryon Density 

The extreme sensitivity of the BBN-inferred estimates of AN, to the adopted helium abundance (and 
its large errors), is responsible for the relatively large error in the BBN-inferred value of AN,. An 
alternate approach avoiding 4 He has been suggested by Nollett and Holder (2011) [I] (see, also, Pettini 
& Cooke (2012) [57j). 

In the best of all worlds the BBN-inferred parameter values should be compared with those in- 
ferred, independently, from the CMB, complemented when necessary to break degeneracies among the 
parameters by other astrophysical data from, e.g., large scale structure, supernovae, and the Hubble 
constant. In the presence of possible new physics, this would enable a probe of the constancy (or 
not) of these parameters in the early Universe epochs from BBN until recombination. However, if 
it is assumed that AN, and 7710 are unchanged from BBN to recombination, the information pro- 
vided by BBN using the helium abundance may be replaced with that from the CMB-determined 
baryon density: 7/io(CMB) = 6.190 ± 0.115 [36]. Using this value in combination with deuterium, 
f?D = Vw — 6(5 — 1) = 5.96 ± 0.28 (for £ = 0), leads to a smaller estimate of AN, but, with a 
larger uncertainty resulting from the much weaker dependence of t/d on AN,: AN, = 0.48^Qg3 
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(N e// = 3.53i^J The 68% and 95% contours for these results are shown in Fig.g While this 
provides some support, once again, for the presence of one sterile neutrino, the absence of dark radi- 
ation (ANj, = 0) is consistent with these results at the ~ 68% confidence level. For these values of 
rjD and 7710 the BBN-predicted helium abundance is Yp = 0.2541 ± 0.0081, very close to(well within 
the errors of) the observationally-inferred value adopted here. Once again, the BBN-predicted lithium 
abundance is a problem: A(Li) = 2.69 ± 0.05. 




Figure 9: The 68% (solid) and 95% (dashed) contours in the A.N U - 7710 plane derived from BBN 
deuterium and the CMB constraint on the baryon density |36j . 

This approach, replacing 4 He with the CMB determined baryon density parameter, could also 
be used to constrain a lepton asymmetry. The corresponding constraint on the neutrino degeneracy, 
£ = — 0.18±0.24, while entirely consistent with £ = 0, has an uninterestingly large uncertainty resulting 
from the very weak dependence of tjd on £. 

0.5.3 Supplementing BBN With The CMB To Constrain £ ^ And AN„ ^ 

The BBN-predicted primordial light element abundances depend on all three of the key parameters 
{7710, ANj,,£}. However, the uncertainty in the observationally-inferred relic abundance of 3 He, along 

For the single, most precise deuterium abundance found by Pettini & Cooke [27|, Vo = 6.10±0.24. If this abundance 
is identified with the primordial deuterium abundance, AN„ = 0.18 ± 0.55 (N e // = 3.22 ± 0.55). 
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Figure 10: The ±lcr band (red) in the £ — AN,, plane from BBN using the D and 4 He constraints. 
The blue band is the ±1 a range for AN^ from Joudaki (2012) [3D] . 



with the lithium problem(s), leaves only two, relatively well constrained primordial abundances, those 
for D and 4 He. From BBN alone and these abundances, all three parameters can't be determined 
independently but, one of them can be eliminated resulting in a relation (degeneracy) between the 
remaining two. For example, t?d and r^He may be used to eliminate 7710, leading to £ = £(AN^; ynp, Yp), 
where 

145£ = 106(5 - 1) + t?d - me- (33) 



This constraint on £ versus ANj, is shown by the red band in Fig. 10 Without an independent 



constraint on ANj,, the degeneracy between £ and AN^ seen in Fig. 10 cannot be broken. The CMB 
provides such a constraint. If BBN, using the observationally- inferred D and 4 He abundances, is now 
supplemented with an independent constraint on AN„ from the CMB, then a combined constraint, 
allowing for both dark radiation and lepton asymmetry, may be found. Using the CMB results from 
WMAP7 [36] . ACT [37] . and the SPT [381 139] , along with complementary constraints from large scale 
structure, the Hubble constant, supernovae, and galaxy clusters, Joudaki (2012) [30] finds ^ e ff = 
3.87 ±0.42, corresponding to S = 1.065 ±0.032. Using this result in Eq.[33] leads to £ = 0.009 ±0.035, 
entirely consistent with £ = 0. The 2a upper bound to |£| is 0.079, corresponding to an upper bound 
to AN„(£) = 0.008. Once again, the BBN-predicted lithium abundance is high, A(Li) = 2.70 ± 0.06, 
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reinforcing the lithium problem. 

A similar approach may be used to eliminate AN^ instead of 7710 to find £ = £(7/10; yDP>Yp) and, 
to use the CMB for a constraint on i]iq. But this approach, which is also consistent with £ = 0, leads 
to a less precise constraint: £ = —0.012 ± 0.052. 

0.6 Summary And Anticipation Of Future Results 

The BBN results presented here provide modest support for the presence of dark radiation (ANj, > 
at ~ 1.5 a), to go along with the more robust evidence for dark matter and dark energy. The current 
observational data (D and 4 He), while allowing for the presence of a sterile neutrino (ANj, = 1 at 
< 1 a), disfavors the presence of two sterile neutrinos (ANj, = 2 at ~ 2.8 a). If, instead, it is assumed 
that AN„ = 0, a small, but non-zero lepton asymmetry is favored, also at the ~ 1.5 a confidence level. 
In contrast, if both AN y and £ are allowed to vary freely and, if BBN (D and 4 He) is supplemented 
by a CMB constraint on AN U , a vanishing lepton asymmetry (£ = 0.009 ± 0.035) is favored. 

Currently, as may be seen from Fig.[7J there is very good agreement between the BBN and CMB 
constraints on the baryon density and dark radiation (the CMB is insensitive to a small or mod- 
est lepton asymmetry). For many years BBN provided the best constraints on the baryon density 
(7/10 or fifth 2 ) and on dark radiation (ANj, or S), as well as the only constraint on lepton asymmetry. 
With WMAP7 [H] and other CMB datasets [371 ESI ES] the best constraints on the baryon density 
now are from the CMB, which allow for a factor of ~ 2 — 3 more precise determination of 7/10 . However, 
in the present, pre-PLANCK era, BBN still provides the best dark radiation constraint, albeit with an 
uncertainty smaller than that from the CMB by only a factor of ~ 1.5 — 2. It is expected that with the 
publication of the PLANCK data the dark radiation torch will pass to the CMB. Depending on what 
PLANCK finds, it may be possible to establish the presence of dark radiation (a sterile neutrino?) at 
the ~ 5cj level if, for example, PLANCK should find, AN^ = 1 ± 0.2. If, however, PLANCK should 
find (the best or worst of all worlds?) AN^ = 0.5 ± 0.2, the presence dark radiation will be favored 
but, that of a sterile neutrino will be somewhat disfavored^ 

For all the possibilities considered here (AN y = £ = 0; £ = 0, AN^ / 0; AN U = 0, £ ^ 0; AN,, / 0, 
£ / 0), the BBN-predicted lithium abundance hardly changed at all (2.68±0.06 < A(Li) < 2.70±0.06). 
This insensitivity is easy to understand since the BBN-predicted lithium and deuterium abundances 
are strongly correlated 

r? Li = r/ D + 3[(5-l)-£]. (34) 

Because the values of 3(S — 1) and 3£ are almost always small compared to t?d, the corrections to a 
perfect lithium - deuterium correlation are generally at only the few percent level. A solution to the 
lithium problem is not to be found with dark radiation or a lepton asymmetry. 

0.6.1 Anticipating The Future 

The future for the key parameters related to the baryon abundance (t/b) and the presence, or not, 
of dark radiation (AN;,) lies with the CMB and the anticipated results from the PLANCK mission. 
While it is impossible to predict the central values PLANCK will find for 7/b or AN^, it is possible to 
forecast the precision to be expected from the PLANCK data and analyses |41l I42j . Such forecasts 
suggest that the uncertainty in the baryon abundance determination will be of order 17(7710) ~ 0.03, 

6 It is perhaps worth recalling that the contribution of a light, thermalized scalar corresponds to AN„ = 4/7 = 0.57. 
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nearly an order of magnitude better than the current BBN precision. The same forecasts suggest that 
AN„ will be constrained to dfANj,) ~ 0.2, or better. This would result in an improvement over the 
current BBN precision by a factor of ~ 2.5. 

If it is assumed that the PLANCK values of the key parameters are identical to those at BBN, 
ignoring their possible evolution from the epoch of BBN until recombination, then these values may 
be used in combination with BBN to predict the relic abundances^] For example, for deuterium, it is 
anticipated that PLANCK will constrain ijd with a precision of <t(??d) ~ 0.1 or, to < 2% for r/D ~ 6. 
The largest uncertainty in the BBN-predicted deuterium abundance at present and in this anticipated 
future arises from uncertain nuclear reaction rates [4J. It can be hoped that this uncertainty may be 
reduced by new laboratory data, reducing the error in the BBN-predicted value of ?/dp by perhaps a 
factor of ~ 2. This would lead to a reduction in the error in the inferred value of t/d by nearly a factor 
of two, <t(t7d) « 0.3 -¥ 0.15. 

For helium, PLANCK may constrain 7/h c to cr(r/Ho) ~ 1-6 which, while still large, is a factor ~ 2.3 
smaller than the current BBN uncertainties. This suggests that using the CMB determined values of 
7710 and AN,,, the BBN-predicted primordial helium mass fraction will be known to c(Yp) < 0.003, 
a precision anticipated to also be attainable in an independent determination of Yp from the CMB 

SB 021- 

For lithium, PLANCK may constrain r/Li to c(?/Li) ~ 0.06 or, to better than ~ 1% for r/Li ~ 6. 
However, as for deuterium, the precision of the BBN-predicted primordial lithium abundance is limited 
by the nuclear physics uncertainties (~ 10%). Nonetheless, it will be very interesting to see if the 
PLANCK data support or, possibly eliminate, the lithium problem(s). 

Although the CMB is insensitive to a lepton asymmetry, as may be seen from Eq. 33 , a combination 
of BBN and CMB constraints on AN^, dbp, and Yp can constrain a neutrino degeneracy, provided 
that the lepton asymmetry is very large compared to the baryon asymmetry. For example, for the 
anticipated CMB constraints on ^(ANj,) ~ 0.2 and on the primordial helium abundance, c(Yp) ~ 
0.003, along with a BBN constraint on y DP , <r s « 0.018 or, a VL » (vr 2 /4C(3))o- ? « 0.036 « 6 x 10 7 r? B - 



0.6.2 Summary 

This review finds itself on the cusp of potentially great changes. Current BBN and CMB data provide 
strong support for the presence of (at least) three SM neutrinos, thermally populated during the 
early evolution of the Universe. This provides indirect support for the so far invisible, relic neutrino 
background. The new CMB and large scale structure data have the potential to constrain the baryon 
asymmetry and the presence, or not, of dark radiation to new levels of precision, testing BBN and the 
current estimates of the relic abundances of the light elements. It will be of great interest to compare 
and contrast the current BBN results with those from the new data and to see what we may learn 
about new physics, including neutrino physics, beyond the standard models of particle physics and 
cosmology. 
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